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^ ' 1 INTRODUCTION 



ABSTRACT 

We present a detailed numerical study of the impact that cosmological models featuring 
a direct interaction between the Dark Energy component that drives the accelerated expansion 
of the Universe and Cold Dark Matter can have on the linear and nonlinear stages of structure 
formation. By means of a series of collisionless N-body simulations we study the influence 
that each of the different effects characterizing these cosmological models - which include 
among others a fifth force, a time variation of particle masses, and a velocity-dependent accel- 
eration - separately have on the growth of density perturbations and on a series of observable 
quantities related to linear and nonlinear cosmic structures, as the matter power spectrum, 
the gravitational bias between baryons and Cold Dark Matter, the halo mass function and the 
halo density profiles. We perform our analysis applying and comparing different numerical 
approaches previously adopted in t he literature, and we address the partial discrepancies re- 
cently claimed in a similar study by Li & Barrow) (1201 Obh with respect to the first outcomes of 
lBaldietalJ ( l2010h . which are found to be related to the specific numerical a pproach adopted 
in the former work. Our results fully confirm the conclusions of lBaldi et al.l (12010 ) and show 
that when linear and nonlinear effects of the interaction between Dark Energy and Cold Dark 
Matter are properly disentangled, the velocity-dependent acceleration is the leading effect 
acting at nonlinear scales, and in particular is the most important mechanism in lowering the 
concentration of Cold Dark Matter halos. 
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The two pillars on which the presently accepted standard cosmo- 
logical model is based, as witnessed by its acronym ACDM, are 
a cosmological constant A which determines an acceleration of 
the expansion of the Universe, and a new type of non relativistic 
massive particles - which go under the name of Cold Dark Mat- 
ter (CDM) - that source the gravitational potential wells in which 
cosmic structures can form. Far from being the ultimate descrip- 
tion of our Universe, the ACDM model represents an efficient way 
to parametrize our ignorance about the fundamental constituents of 
roughly 95% of its energy content. 

Although the combination of a cosmological constant A 
and of a significant fraction of CDM particles provides a very 
goo d fit to most o f the prese n tly available observations (see 
e.g. iKoma tsu et aL 2009, '20 lO'; 'Percival et al.' '200 1"; 'Cole et alJ 



2005 ; Reid et al.. .2010.: .Riess et al, .1998.: Perlmutter et al.. ,199' 



Astier et al.ll2006l : iKowalski et al.ll2008riPercival et al.ll2010l) . the 

fundamental nature of the two main components of the Universe 
remains an open question. On one side, in fact, the observed value 
of the energy scale associated with the cosmological constant dif- 



fers by tens of orders of magnitude from its theoretical predictions, 
thereby making of the cosmological constant an extremely fine- 
tuned parameter of the model. On the other side, all the proposed 
CDM candidates have so far evaded any attempt of direct detection 
and their properties can still be inferred only by cosmological and 
astrophysical observations. 

It is in this context that alternative models have been pro- 
posed, by taking the standard ACDM as an asymptotic state of 
more complex underlying scenarios. These generally involve, as 
a first step, the promotion of the cosmological constant A to a 
dynamical quantity - dubbed dark energy (DE) - which evolves 
during the expansion history of the Universe just as all its other 
physical constituents. A convenient way to represent this sce- 
nario involves the dynamical evolution of a classical scalar field 
in a self interaction potential, which goes under the name of 
quintessence ( Wetterich 1988; Ratra & Peebles|[T988h or k-essence 
dArmendariz-Picon et al.ll2001l) . 

Other possible alternatives with respect to the standard 
cosmological scenario are given by modifications of the laws 
of gravity at large scales or at low spatial curvatures (e.g. 
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IStarobinskvl [l980l ; IHu & Sawickil l2007h . or by relaxing the as- 
sumpdons of homogeneity and isot ropy of the Universe (see e.g. 
iGarcia- Bellido & Haugboelldl2008l) . 

Particular attention has been devoted, in recent years, to the 
idea that a dynamical DE scalar field might have direct interac- 
tions (besides gravity) with other cosmic fluids, by directly ex- 
changing energy- mom entum (Wetterich 1995; Amendola 2000; 
lFarrar & Peebles! I2OO4I ; Umendola et alJi2008D . While such direct 
interactions with standard model particles would be heavily re- 
stricted by available observational constraints, the same does not 
happen for the case of a selective interaction with CDM particles 
only, as first suggested by Damour et al. ( 1990), for which obser- 
vational bounds are much weaker. 

These models have been extensively studied in the recent past, 
with a particular focus on their possible distinctive effects on struc- 
ture formation. In fact, all interacting DE models predict the exis- 
tence of an additional attractive force - mediated by the DE scalar 
field - between massive particles which is expected to alter the 
growth of density perturbations. Besides this "fifth-force", cou- 
pled massive particles also experience, in interacting DE models, 
a "modified inertia" determined by the energy-momentum transfer 
between the DE and the CDM sectors. 

The combination of these new physical effects and of the 
modified cosmic expansion - that arises as a consequence of 
the dynamical nature of the DE field - determines a modi- 
fied evolution of density perturbations thereby providing possi- 
ble observational signatures of the models. The effects of inter- 
acting DE models on the evolution of cosmic structures have 
been stud ied both m the Imear regime (Amendola 2000, 2004'; 



where the interaction with CDM particles is represented by a source 
term in the respective continuity equations for the two fluidj^: 



Pettorino & Baccigalupi 2008; Di Porto & Amendola 2008; Baldi 
2010! ) and in the nonlinear regime by relying on simplified configu- 



rations, as for the case of spherical c ollapse tMainini & Bonomettol 
l2006l ; IWintergerst & Pettorinol2010l). or in ful l gener a lity by means 
of N -body simulations (e.g. by iBaldi et alj 1 20 id ; iBaldil I2OI0I ; 
iLi & Barrow 2010a). 

In this work we present an extension of previous analyses 
about how the different new physical effects featured by interact- 
ing DE models individually contribute to the overall modifications 
of the properties of linear and nonlinear cosmic structures, and we 
try to clarify some confusion that has been made in the literature 
concerning the relative importance of the different effects in these 
two different regimes. 

The paper is organized as follows: in Sec.[2]we briefly review 
the main equations of interacting DE models both concerning the 
background evolution and the growth of linear perturbations, with 
a particular focus on the issue of model normalization which will 
be central for the rest of the present study; in Sec. [3] we discuss 
the different effects that contribute to a modification of newtonian 
dynamics at small scales for CDM particles; in Sec.|4]we present 
our set of simulations, describing which effects have been included 
in each simulation, and how the relevant quantities have been nor- 
malized in order to allow a meaningful comparison of the differ- 
ent outcomes; in Sec. [5] we present our results concerning linear 
and nonlinear properties of cosmic structures under different im- 
plementations of the interacting DE effects; finally, in Sec. |6]we 
draw our conclusions. 



-P{<f>)<f>Pc 
+/3(<^)#c 



(1) 

(2) 



where pc is the CDM density, p^ = <^^/2 + Vijj)) is the DE den- 
sity, H is the Hubble function, and where an overdot represents a 
derivative with respect to cosmic time t. The function deter- 
mines the strength of the DE-CDM interaction and together with 
the scalar potential V[cj>) fully specifies the model. 

Interacting DE models described by Eqs. ( Ill2t . where all 
the other cosmic fluids (as radiation, baryonic matter, mas- 
sive neutrinos, etc..) obey standard continuity equations, have 
been widely studied in the p ast both f or the case of a con- 
stant coupling (as e.g. by WetterIci][T993 ; lAmendolall2000l. |2004| ; 



IPettorino & BaccigalupilbOOsl . and references therei n) and for th e 
more general situation of a field depentend coupling l lBaldill2010h . 
In the present work we will focus on the former situation, and 
we will therefore assume P{<f>) = /3 for the rest of the paper. 
We will also always assume an exponential form for the potential 
1/(0) oce"""^, with a = 0.1. 

The background evolution of constant coupling models is 
characterized by a scaling regime during matter domination where 
the two interacting fluids (DE and CDM in our case) share a con- 
stant ratio of the total energy budget of the universfl (known as 
the (jiMDE regime), thereby determining the presence of a sizeable 
fraction of Early DE (EDE hereafter) which alters the standard ex- 
pansion history with respect to a corresponding ACDM model with 
the same set of cosmolgical parameters at the present time. This 
scaling regime is sustained (for the case of positive couplings on 
which we focus in the present work) by the energy transfer from 
the CDM particles to the DE scalar field, which determines in turn 
a decrease in time of the mass of CDM particles, according to the 
modified continuity equation (TJ. Therefore, a modified expansion 
rate at z > and a time evolution of CDM particle masses are two 
common features of any interacting DE model which can affect the 
growth rate of density perturbations. 

For the latter effect, it is clearly important the normalization 
of the CDM masses: if one wants to compare models that share the 
same cosmological parameters at 2 = 0, then the mass of CDM 
particles within interacting DE models will be larger in the past, 
corresponding to an effectively higher value of the CDM density 
Pc- However, it is important to stress here the distinction between 
the mass variation, i.e. the fact that rhc 7^ 0, and the mass nor- 
malization, i.e. the effective pc{z) during the expansion history of 
the universe. This distinction is particularly important for what we 
will show in the rest of the paper, since the former effect is found 
to have significant implications for the nonlinear regime of struc- 
ture formation and for the internal dynamics of collapsed objects, 
while the latter is primarily affecting the linear evolution of den- 
sity perturbations. Not recognizing this distinction can cause some 
confusion in the study of the nonlinear effects of interacting DE 
models, as we will discuss below. 

Besides these two features, the stud y of the linear perturba- 
tions evolution in inte racting DE models (Amendola '200dl. l2004l; 
IPettorino & Baccigalupi 2008) has led to the identification of other 
two relevant effects which directly influence the processes of struc- 



2 INTERACTING SCALAR FIELDS 

We consider interacting DE models where the role of DE is played 
by a scalar field (f) evolving in a self-interaction potential V((l)), and 



^ Throughout the paper we will use units in which the reduced Planck mass 

is assumed to be unity, Mpi = I/VSttG = 1. 

^ We always assume spatial flatness such that ptot = 3H^. 
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ture evolution. The dynamic equation for CDM density perturba- 
tions in interacting DE scenarios 

5^+(2H ^ l3,},'^S,^^H^{l + 2l3'')njc + ntSb] =0, (3) 

shows in fact also the presence of an additional friction term di- 
rectly proportional to the coupling 

- P^S, , (4) 

and of an effective enhancement of the gravitational pull for CDM 
fluctuations by a factor (1 + 2/?^) which is known as the "fifth- 
force". Both the extra friction term and the fifth-force accelerate 
the growth of CDM density perturbations in the linear regime as 
clearly shown by Eqn. 

2.1 Model normalization 

It is of particular relevance for the kind of analysis carried out in the 
present work to clarify which normalization is assumed for the dif- 
ferent cosmological models under investigation (that in this case are 
a ACDM cosmology and an interacting DE scenario with constant 
coupling /3) and for the different numerical realizations of the latter 
in the test simulations that we are going to use for our discussion. 
The normalization of the models concerns both the background and 
the linear perturbations evolution in the different test cases, which 
we separately discuss here. 

2.1.1 Background normalization 

In the present work we adopt the same type of nor malizat i on as- 
sumed in the previous works o f lBaldietaUhoiOh and lBaldiUloToh 
for the background evolution, where both the ACDM and the inter- 
acting DE model share the same background cosmological param- 
eters at 2: = 0. This is realized by assuming a specific set of cosmo- 
logical parameters (given by Ho, ^Ic, ^t, ^r, ^de) in accordance 
with the latest observ ational constraints from the WMAP satellite 
jKomatsu et"^l2010 l) and by then integrating backwards in time 
the full system of background equations for each specific model. 
As a consequence of this procedure, both the ACDM and the inter- 
acting DE model will have exactly the same cosmological parame- 
ters at 2 = and will therefore differ from each other at higher red- 
shifts. In particular, the mass of CDM particles will be larger at high 
redshifts due to the interaction with the DE scalar filed (j) and to the 
continuous flow of energy from the former to the latter. Due to the 
existence of the 0MDE scaling regime, where a sizeable fraction of 
EDE is present during matter domination, and to the requirement of 
spatial flatness (Jlc+fib+f^r+fiDE = 1), if a larger fraction of DE 
is present in interacting DE models during matter domination, then 
the total value of f2c will necessarily have to be smaller at the same 
time as compared to ACDM. A larger value of pc and a smaller 
value of Sic = pc/SH^ then obviously implies a larger value of the 
Hubble function during (/>MDE with respect to ACDM. This is what 
is shown in Fig. [T] where the evolution of the cosmological den- 
sity parameters and of the Hubble function are shown for ACDM 
and for our interacting DE model. This is consistent with most 
of the other studies on interacting DE and EDE models (as e.g. 
Doran et al. 2001; Wetterich 2004: Pettor ino & BaccigalupilbOOSl : 
Francis etal. 2008 : Grossi & Springel 200^ . 

This type of normalization therefore allows to compare differ- 
ent models that have exactly the same cosmological parameters at 
z = 0, thereby providing a homogeneous set of models to be con- 
fronted with each other. Other studies have adopted different types 



of normalization (as e.g. in iLi & Barrowll2010al lbh where different 
interacting DE models are normalized by assuming the same ini- 
tial conditions in the early universe for the background dynamical 
quantities and by evolving these initial conditions forward in time 
with a trial and error procedure until suitable solutions are found. 
This procedure, however, clearly does not guarantee to get the same 
values for the background cosmological parameters at z = 0. 

2.1.2 Linear perturbations normalization 

A completely independent choice of normalization between differ- 
ent interacting DE models concerns the amplitude of linear den- 
sity perturbations. As a consequence of the new physical pro- 
cesses and of the different background evolution, in fact, inter- 
acting DE models will have a significantly different growth rate 
of density perturbations as compared to ACDM, whe re the devi- 
ation will depend on the strength of the coupling (see lBaldi et al.l 
l2010l : lDi Porto & Amendoiall2008l , for a detailed description of the 
growth factor in interacting DE scenarios). When setting up the ini- 
tial conditions for N-body simulations, one is therefore allowed to 
choose whether to set the desired value of the linear density per- 
turbations at high redshifts, by fixing the amplitude according to 
the observed value of the scalar amplitude .4s from CMB data, 
or at the present time, by assuming the present value of as as de- 
termined by low-redshift measurements as normalizations. In the 
former case, different models will start with identical density am- 
plitudes at early times a nd will en d up with different values of as at 
z = (as shown e.g. in lSaldi & Pettorino 2010), while in the lat- 
ter the amplitude of the initial density fluctuations will have to be 
scaled with the appropriate growth factor of each speci fic model in 
order to give the same value of ag at pre sent (as done in iBaldi et al] 
120101 : lBaldill2010l : lLi & Barrowll2010ah . 

Both types of normalization have been used in the past in dif- 
ferent numerical studies aimed at different types of analysis of non- 
linear structure formation within interacting DE models. However, 
the normalization at z = is clearly the preferred choice if one 
aims at directly comparing the effects of interacting DE on the non- 
linear dynamics of CDM particles within collapsed structures, since 
with such normalization all the models will share at the present time 
the same cosmological parameters and the same linear density per- 
turbations normalization, and all the differences in the properties 
of low redshift structures will be related to the nonlinear behavior 
only. 

It is nevertheless important to stress here that if any of the new 
physical effects of interacting DE is artificially suppressed in an N- 
body run, then the effective linear growth factor in the simulation 
will change and a different value of erg at z = will be obtained 
even if the initial conditions were normalized with the preferred 
procedure described above. This might provide interesting infor- 
mation on the impact that the specific effect which has been sup- 
pressed has on the linear regime of structure formation, but could 
also introduce some confusion in the study of the nonlinear regime, 
since nonlinear structures forming in such test simulation will be 
embedded in a linear density field with a different normalization. 



3 MODIFIED GRAVITATIONAL DYNAMICS 

If we want to follow the evolution of density perturbations beyond 
the linear regime and be able to predict the features that interacting 
DE imprints on the highly nonlinear objects that we can directly 
observe in the sky, Eqn. ^ is no longer sufficient and we need to 
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Figure 1. (Color online) Left panel: The background evolution of a standard ACDM model and of an interacting DE model with j3 = 0.24. In this example 
we have not included the uncoupled fraction of baryonic matter to simplify the plot. The figure shows the evolution of the dimensionless density in radiation 
(green), matter (black) and dark energy (red) as a function of the e-folding time (defined as the logarithm of the scale factor a) for ACDM (solid) and interacting 
DE (dashed). The smaller plot shows the ratio of the same quantities to the ACDM case. Right panel: The ratio of the Hubble function of the same two models 
with respect to ACDM as a function of the e-folding time. A larger Hubble function with respect to ACDM is consistent with all previous studies on interacting 
DE and EDE, as explained in the text. 



rely on numerical integrations. In order to do so, one has to identify 
how the interaction between DE and CDM affects the laws of new- 
tonian dynamics that govern the evolution of structure formation 
in the newtonian limit of General Relativity and implement these 
effects in to N-body a lgorithms. 

Baldi et alj ( 12010.) have shown that the acceleration equation 
for a CDM particle in an interacting DE cosmology for the case 
of a light scalar field (i.e. a scalar fi eld model for wh ich = 
d^V/A(f)'^ < H) takes the form (see lBaldi et ai]|2010l for the de- 
tails of the derivation): 



G(l + 2/32)mjri. 



(5) 



where Vi is the velocity of the i-th particle, r^j is the vector dis- 
tance between the i-th and the j-th particles, and the sum extends 
to all the CDM particles in the universe. 

By having a look at Eqn. ([5}, one can clearly identify the same 
coupling-dependent terms already encountered for the linear per- 
turbations equation ([S}. The friction term of Eq. (O now appears as 
a "velocity-dependent" acceleration 



a„ = I3<j)v 



(6) 



which depends on the vectorial velocity of CDM particles, while 
the "fifth-force" term appears in the same form as in Eq. l|3j, i.e. as 
an effective rescaling of the gravitational constant G between CDM 
particles pairs by a factor (1 + 2/3^). 

It is however very important to notice a striking difference be- 
tween the linear and the nonlinear regimes, that will have signifi- 
cant implications for the discussion carried out in this paper: while 
in the linear regime the friction term always accelerates the growth 
of structures, in the nonlinear case this depends on the relative ori- 
entations of the velocity and of the gravitational acceleration of 
each CDM particle. As a consequence, a particle moving towards 
the local potential minimum will experience an effectively larger 
potential gradient while a particle moving away from the local po- 
tential minimum will conversely feel an effectively smaller poten- 
tial gradient. For the realistic situation of nonlinear virialized ob- 
jects, where tangential velocities are non negligible with respect to 



radial velocities, the velocity-dependent acceleration will therefore 
have a completely different effect than in the linear regime. 

For this reason, when comparing the properties of nonlinear 
structures one should avoid to consider the linear friction term 
ad the nonlinear velocity-dependent acceleration as a single phe- 
nomenon, and always distinguish between its linear and nonlinear 
behavior. As we will show below, not making this clear distinction 
between the two regimes can cause some further confusion in the 
determination of which are the most relevant effects of interacting 
DE for the nonlinear dynamics of CDM particles. 



4 SIMULATIONS 

The study of the nonlinear effects of interacting DE models has 
been approached in recent years by means of suitably modified N- 
body algorithms. The first hydrodynamical high-resolution N-body 
simulations of interacting DE models have been perform ed with 
a mo dification of the parallel TreePM c ode GADGET-2 ( iSpringell 
l2005l) and presented in lBaldi et al.l j20I(]|) . Other studies have then 
been carried out by means of mesh or Tree based N-body algo- 
rithms, but without h ydrodynamics (see e.g. ,Li & Ba rrow 2010^; 
iHellwing et"ai]|2010l) . and showed a good agreement with the first 
results of lBaldi et all j2010l) . 

These could be summarized as follows: 



• The DE-CDM interaction determines a faster growth of lin- 
ear density perturbations as compared to ACDM, such that in order 
to obtain the same normalization of the power spectrum amplitude 
at z = (i.e. models with the same ag, at the present time) it is 
necessary to rescale the density fluctuations amplitude in the ini- 
tial conditions with the appropriate growth factor for each different 
interacting DE model; 

• The selective interaction only between DE and CDM leaves 
the baryons completely uncoupled, and as a consequence baryonic 
and CDM density fluctuations evolve with a different rate in in- 
teracting DE models; this determines a significant reduction of the 
halo baryon fraction of collapsed objects at z = 0; 
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• If models are normalized at 2: = for what concerns the am- 
plitude of their linear power spectrum, the mass function of halos 
has a comparable shape and amplitude at 2 = for all the models, 
while it differs at larger redshifts for different values of the cou- 
pling /3; this does not hold if the models are normalized at high 
redshift (as e.g. with a normalization at CMB), in which case inter- 
acting DE models show a significant increase of massive ob jects at 
any r edshift as compared to ACDM, as shown bv lSaldi & Pettorinol 

mlo.) : 

• For the case of constant couplings considered in this work, the 
CDM density profiles of massive halos at 2 = are always less 
concentrated in interacting DE scenarios as compared to ACDM; 
this does not necessarily hold for th e more general case of time 
dependent couplings (see lBaldj||2010h . 

Due to the complex interplay of all the new physical effects 
arising in interacting DE models and discussed in Sec.|2]and[3l it 
is an interesting exercise to try to isolate the individual impact that 
each of these different modifications of the standard gravitational 
dynamics has on the final evolution of linear and nonlinear struc- 
tures and to investigate how they separately contribute to each of 
the peculiar features of interacting DE models listed above. This 
kind of analysis had already been carried out in iBaldi et al.l JioToh 
for the effect of concentration reduction in interacting DE models, 
by artificially switching off in the N-body integration some of the 
new physical effects during the latest stages of structure formation 
where nonlinear processes take place. This has been done by com- 
puting the average formation redshift 2form of the most massive 
halos identified in the simulations and by selectively suppressing 
each specific effect of the DE-CDM interaction only after a red- 
shift 2 ni somewhat larger than 2form- For instance, in lBaldi et al.l 
the average formation redshift in the sample of the 200 most 
massive halos formed in the different simulations was found to be 
ZioTui ~ 0.8, and the swhitch-off redshift for the test simulations 
was chosen to be 2ni = 1.5. 

This procedure ensures to isolate as much as possible the non- 
linear effects of the interacting DE models from the linear evolution 
of the overall power spectrum amplitude and therefore to be able to 
compare the properties of nonlinear objects that form in similar 
environments and with comparable formation historie s. The out- 
comes of the analysis carried out with this approach in lBaldi et al.l 
( l2010h showed that the reduction of halo concentrations is primar- 
ily determined by the velocity-dependent acceleration {SJ with a 
minor contribution from the effect of mass variation of CDM parti- 
cles. 

A more recent study carried out by iLi & Barrow! ( l2010bh 
has shown some partial discrepancies with the early results of 
iBal di et al. (2010) claiming that the mass variation effect be ac- 
tually the most relevant mechanism to deter mine the reduction of 
halo concentrations. However, the analysis o f lLi& Ban-owl ( l2010bh 
did not adopt the above mentioned procedure aimed at distinguish- 
ing the linear and nonlinear impact of the different effects of the 
DE-CDM interaction, but rather switched off each individual effect 
right from the start of the numerical integrations. This might deter- 
mine some confusion in the interpretation of the results, and drive 
to misleading conclusions, as we will show below. 

In the pre s ent w ork we significantly extend the analyses done 
bv lBaldi et al] ( 1201 Q) and by Li & Barrow. (201 Obi) and we try to 
address the claimed discrepancies between the two treatments, clar- 
ifying the confusion that might arise when the linear and nonlinear 
effects of interacting DE models are mixed with each other. 

To this end, we will make use of the same code presented in 



lBaIdietalJ ( l2O10l) , to which we refer the interested reader for an ex- 
tensive discussion of the implementation strategies. With this code 
we run a series of test simulations for an interacting DE model with 
/3 — 0.24 and with different numerical setups, and for a fiducial 
ACDM modefl Following the approach discussed in Sec. 12.1.11 
all the models are normalized to have the same set of cosmological 
parameters at 2 = in accordance to the most recent results from 
WMAP iKomatsu et al.ll2010l) . A coupling as large as /3 = 0.24 
is an extreme case and is already ruled out by several independent 
observational constraints (see e.g. [ Bean et al.l2008l : lLa Vacca et al.l 
l2009l : IXij2009l : lBaldi & Viell2010l) . but our aim here is to amplify 
the impact of the interaction in order disentangle more clearly the 
contributions arising from the different specific effects of interact- 
ing DE. This means that the simulations presented in this work are 
not intended to represent realistic models of the universe, but rather 
to provide a clear description of how interacting DE affects struc- 
ture formation in the linear and in the nonlinear regimes. 

Our simulations consist of a cosmological box of 80 comov- 
ing Mpc/h aside filled with 256^ particles for CDM and baryons 
for a total of = Nc + Nf, x 3.4 x 10^ simulation particles. 
This setup determines a spatial resol ution of » 7 k pc/h, which 
is lower with respect to the work of iBaldi et alj ( 12010,) but com- 
parable with iLi & Barrow! ( !2010b!) . The baryons are included to 
give a correct representation of the effective coupling, but in or- 
der to avoid spurious effects besides the ones determined by the 
DE-CDM interaction we do not include hydrodynamical forces in 
the simulations. In other words, the baryonic particles are treated 
as an additional family of uncoupled collisionless particles which 
contribute to reduce the total effective coupling with respect to the 
unrealistic case where the DE couples to all massive particles in the 
universe. The mass resolution is nic (2 = 0) f« 2.0 X 10'-* Mo/h 
for the CDM particles and mj, « 3.9 x 10* Mq /h for the baryons. 

With these numerical parameters we run a series of 8 sim- 
ulations in which we artificially switch off one of the characteri- 
stic effects of inte racting DE, adopting both the procedures used by 
!Baldi et alj bOlOl) and bv !Li & Barrowl ( !2010j) in order to compare 
the outcomes of these test runs with the full interacting DE model 
and with the fiducial ACDM cosmology. It is of particular relevance 
for our discussion to notice here that the former approach ensures 
to only slightly perturb the normalization of the background cos- 
mological parameters and of the linear perturbations amplitude of 
the different test runs with respect to the fiducial ACDM scenario at 
2 = 0. On the contrary, the latter procedure significantly changes 
the normalization of the background cosmological parameters at 
2 = 0, and more importantly shifts by a large offset the linear am- 
plitude of density perturbations in the different runs, such that each 
of the test simulations will have a significantly different value of as 
at 2 = 0, thereby making quite difficult a direct comparison of the 
nonlinear properties of the collapsed objects forming in the differ- 
ent runs. The details of these simulations are described in Table|4l 

Simulations S2- S5 correspond to the procedure adopted by 
!Li & Barrowl ( !2010d) . to which we add simulation S6 where the 
mass variation of CDM particles is also suppressed during the 
whole run but the mass is normalized to the value it has at 2 = 
in the other simulations. These two different normalizations of the 



Please notice the convention assumed for the definition of the coupling, 
which differs from what assumed in part of the literature. A coupling of 
j3 = 0.24 with th is convention corresponds to /3 = 0.3 in the notation of 
lAmendoljJioOOl) . 
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Name 


iVIOClcl 


Modified Hubble function 


fifth-force 


velocity-dependent acceleration 


mass variation 


mass normalization 


so 


ACDM 


inactive 


inactive 


inactive 


inactive 


•mc(z) = m^{Qi)\/z 


SI 


P = 0.24 


active 


active 


active 


active 


mc{z) > mc(0)Vz > 


S2 




inactive 


active 


active 


active 


mc(z) > mc(0)yz > 


S3 




active 


inactive 


active 


active 


mc{z) > mc{0)Vz > 


S4 




active 


active 


always inactive 


active 


mc{z) > mc(0)Vz > 


S5 




active 


active 


active 


always inactive 


mc{z) = mc{oo)\/z 


S6 




active 


active 


active 


always inactive 


mc{z) = rricioyiz 


S7 




active 


inactive for ^ < 2 


active 


active 


mc{z) > mc{0)Vz > 


S8 




active 


active 


inactive for z < 2 


active 


mc{z) > mc(0)V2 > 


S9 




active 


active 


active 


inactive for z < 2 


(mdz) > mc{2)Vz > 2 
\ mdz) = mc{2)Vz <= 2 



Table 1. Table of the different simulations presented in this work. SO cor responds to the standard ACDM model, S 1 to the full interacting DE simulation. 
Simulations S2-S6 adopt the same procedure used by .Li & Barrow! i2010bl) while simulations S7-S9 follow the method o f iBaldi et aljj2010l) . 
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Figure 2. (Color online) The evolution of CDM particle mass as a function 
of redshift in all the simulations of Table |4] The curved line correponds to 
the evolution in the full interacting DE model. This evolution can be stopped 
at different times in the different simulations according to the different nu- 
merical methods discussed in the text. 



CDM particle mass will determine a huge difference in the final re- 
sults, as we will show below. In fact, simply switching off the mass 
variation of CDM particles from the starting redshift of the simula- 
tions but leaving the particle mass normalization to the same initial 
value of the full simulations would result in a much larger value of 
f2c at 2 = as compared to the fiducial ACDM model, while this 
would not happen if the mass is normalized to the value it takes 
at 2 = in the other runs. It therefore comes as no surprise that 
in the former case the global amplitude of linear density perturba- 
tions will grow significantly faster than in the full interacting DE 
model represented by the simulation SI, and significantly slower in 
the latter one. The situation is well explained in Fig. |2l where the 
mass evolution of CDM particles is shown for all the simulations 
of Table|4] Clearly, adopting a mass normalization to the value that 
the CDM mass has at 2 = 60 in t he full model (as don e in the 
simulation S5 and in the analysis bv lLi & Barrowl ^201 Obi) ) forces 
a significantly larger Q,c during the whole expansion history of the 
universe, and as a consequence a much faster linear growth of per- 
turbations that will influence the subsequent nonlinear evolution 
of massive halos. For instance, with this approach halos will form 
much earlier and in a universe with a much larger CDM density 
than they would do in the full simulation. It is therefore again not 
su rprising that h alos forming in this kind of simulation, as found 
bv lLi & Barrowl fcoiObh . will have a larger characteristic overden- 



sity and a somewhat larger halo concentration. However, this is just 
an artifact arising from the different effective cosmological param- 
eters that are forced into the simulation by normalizing the mass 
at high redshifts, and not by the direct effect of the mass variation 
onto nonlinear structures. In other words, this effect will be primar- 
ily a consequence of the different mass normalizations rather than 
of the mass constancy in time. 

The same argument applies to the case of simulation S6, where 
the mass is normalized to the value of mc{z — 0) in the full simu- 
lation. In this case structures will evolve with a lower effective Qc 
than in the full simulation, that will determine a slower growth of 
density perturbations. As a consequence, starting from a lower am- 
plitude with respect to ACDM in the initial conditions, the CDM 
density perturbations of simulation S6 will not have the time to 
catch up the same as of the other models al z — 0, thereby deter- 
mining a later formation time of massive halos. 

A similar situation concerns the suppression of the velocity- 
dependent acceleration. As we explained above, this term has com- 
pletely different effects in the linear and nonlinear regimes of struc- 
ture formation. In the former situation, in fact, this term always ac- 
celerates the growth of linear density perturbations thereby acting 
in the same direction of a larger value of Qc, while in the latter 
it can significantly reduce nonlinear overdensities at small scales. 
This means that suppressing this term right from the beginning of 
the simulations has a similar effect as what discussed above for the 
mass normalization, and could lead to a different amplitude of lin- 
ear perturbations at 2 = as compared to the fiducial ACDM and 
the full interacting DE model. 

On the contrary, suppressing the mass variation and the 
velocity-dependent acceleration only at recent epochs (in this anal- 
ysis we have chosen 2„; = 2), where the nonlinear stages of struc- 
ture formation take place (as done for simulations S8 and S9) en- 
sures that the CDM density and the amplitude of linear density per- 
turbations at 2: = be comparable with both the ACDM reference 
model and the full interacting DE sim ulation. This was in fact the 
procedure adopted by the first work of lBaldi et al.) ( l2010h . and al- 
lows a direct comparison of the impact that each different effect 
of the DE-CDM interaction has on the nonlinear regime of struc- 
ture formation only, without spurious effects induced by different 
background or linear evolutions of the models. 

It is nevertheless important to stress here that this type of study 
is only a way to better understand how the DE-CDM interaction 
acts on the linear and nonlinear stages of structure formation, while 
it does not represent any realistic situation. The only realistic sce- 
nario is the full model where all the above mentioned effects are 
simultaneously and consistently included. For instance, the mass 
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variation and the velocity-dependent acceleration are two aspects of 
the same phenomenon, namely the conservation of particles' mo- 
mentum, and none of the two can be in place without the other. It 
is nevertheless interesting to perform this kind of analysis to get a 
deeper insight into the dynamics of coupled systems. 

In the next Section we will present the outcomes of the whole 
set of simulations of Table|4]both concerning the linear and nonlin- 
ear stages of structure formation, and we wil l show that the asse ssed 
partial disagreement between the results of lBaldi et all ilOld) and 
lLi&Barro\^J ( l2010bl) is actually just a consequence of the proce- 
dure adopted in the latter study for the suppression of the individual 
effects under investigation, which inevitably mixes linear and non- 
linear effects of the DE-CDM interaction. 



5 RESULTS 

We now compare the results of our test simulations with the refer- 
ence ACDM model and with the full implementation of interacting 
DE carried out in simulations SO and SI, respectively. We consider 
several different observable quantities related to the linear and the 
nonlinear evolution of structures and we try to identify which of the 
various effects of the DE-CDM interaction have a more prominent 
impact on each of these observables. 



5.1 Matter Power Spectrum 

As a first step, we compute for all of our simulations the total matter 
pawer spectrum P{k, z) as a function of inverse scale k at different 
redshifts z, and the separate power spectra for CDM and baryonic 
particles, Pc{k,z) and Pb{k,z), respectively. In Fig. [3] we show 
the total matter power spectrum P{k, z) at redshifts z = 10, 1, 0. 
The upper panels refer to the simulations where the different in- 
dividual effects of the DE-CDM interaction have been suppressed 
fr om the start of the sim ulation at Zi — 60 (as done in the analysis 
of lLi & Barro"wll2010bl) . while the lower panels refer to the proce- 
dure of switching off the specific effects only at late times, i.e. at 
2 < 2ni = 2, as for the study of iBaldi et al.l feOld) . 

First of all, one can clearly identify how the DE-CDM interac- 
tion affects the total matter power spectrum by comparing the black 
(ACDM) and the red (interacting DE) solid curves in the z = 
plots. While the two power spectra have (by construction) the same 
amplitude at the largest scales, i.e. are normalized to give the same 
linear perturbations amplitude at z = 0, at smaller and smaller 
scales (i.e. larger values of k), the interacting DE model features a 
progressively stronger lack of power as compared to ACDM. This 
already shows how the linear and the nonlinear regimes are very 
differently affected by the DE-CDM interaction, since nonlinear 
scales show a significant reduction of power while linear scales re- 
main practically unaffected. 

We can now compare how the different effects of the DE- 
CDM interaction individually contribute to the total distortion of 
the matter power spectrum by having a look at how suppressing 
each specific effect influences the power spectrum at different z. 
Already at a first glance, it is easy to realize how the top panels 
of Fig.Opresent a much larger scatter of the different curves at all 
scales, including the largest scales available in the simulation box, 
as compared to the bottom panels. This already clearly shows how 
suppressing any of the characteristic effects of the DE-CDM inter- 
action from the starting of the simulations (as done in the upper 
panels) determines a large change in the respective values of erg at 



the present time. In particular, suppressing the background modifi- 
cation of the expansion history (i.e. replacing the proper expansion 
history of the interacting DE model with the standard ACDM one) 
produces a faster growth of density perturbations that consequently 
reach a higher normalization at the present time, as shown by the 
cyan line in the top panels of Fig. [3] This is fully consistent with 
the fact that the interaction determines a larger value of the Hub- 
ble function al z > (as shown in the right panel of Fig.[TJ, and 
therefore replacing it with the ACDM Hubble function determines 
a slower expansion rate as compared to the full simul ation. This re- 
sult is however exactly the opposite of what shown bv Li & Barrowl 
f2010b), where a different normalization has been used and where 
interacting DE models feature a slower expansion rate than ACDM 
at 2 > 0. 

The suppression of mass variation, as we discussed above, also 
determines a different effecive Q,c{z) throughout the simulation, 
thereby inducing a faster or slower growth of density perturbations 
depending on whether the constant mass value is fixed to the initial 
(rric — mc{zi), simulation S5, dotted orange line) or to the final 
(rric = mc(0), simulation S6, dashed orange line) value taken in 
the full simulation SI, respectively. This effect clearly appears in 
the upper panels of Fig. [3] where the high- and the low-redshift 
normalizations of the CDM particle mass give rise to higher and 
lower values, respectively, of the power spectrum amplitude at all 
scales. 

The suppression of the fifth-force (green lines) and of the 
velocity-dependent acceleration (blue lines) produce weaker effects 
on the overall amplitude of the power spectra, while it is already 
clear from the upper panels of Fig[3] that the velocity-dependent 
acceleration has a stronger impact at nonlinear scales since its sup- 
pression (simulation S4) clearly shows an increase of power at large 
values of k with respect to the full simulation S 1 which is not de- 
tected if the fifth-force is suppressed (simulation S3). 

The analysis of the nonlinear effects becomes then much more 
clear if we move to analyze the outcomes of test simulations run 
adopting the procedure of Baldi et al.H2O10l) . As we stressed above, 
this procedure seems more suitable if one wants to compare the 
nonlinear effects of interacting DE models since it does not alter the 
overall linear normalization of the different test simulations under 
study. This is very clear by looking at the lower panels of Fig. [3] At 
z > Zni all the simulations show the same power spectrum at all 
scales, as expected, with a lower overall amplitude as compared to 
ACDM due to the scaling of the initial conditions. At z ^ z^i the 
different test simulations start to deviate from each other at small 
scales, while the large scale normalization stays the same until z = 
0, where it finally coincides also with the ACDM one. 

This allows to compare the importance of the different effects 
of the DE-CDM interaction on the nonlinear scales for structures 
embedded in a density field with the same linear amplitude normal- 
ization at large scales. This more meaningful comparison already 
shows a clear hierarchy of the relative contribution of the different 
effects to the total distortion of the nonlinear matter power spec- 
trum. If the fifth-force (green) seems to have basically no effect 
in this context, the velocity-dependent acceleration (blue) is clearly 
the most important effect in erasing power at small scales, while the 
mass variation (orange) has a sizeable b ut subleading impact . This 
result is in contrast with what shown bv lLi & Barrowl ( l2010bl) , and 
shows how a high-redshift normalization can determine a poten- 
tially misleading superposition of linear and nonlinear effects. 



© 2010 RAS, MNRAS 0Q0.[Tlfl4l 



8 M. Baldi 




Figure 3. (Color online) The matter power spectrum of ACDM (black) and interacting DE (red) as a function of inverse scale k at different redshifts. The upper 
panels show the test simulations run suppressing specific effects of the DE-CDM interaction right from the start of the simulations. Lower panels show the 
results of simulations with a suppression only during the latest stages of structure formation, i.e. for z z^i = 2. Clearly, the large scatter present in the upper 
panels witnesses the superposition of linear and nonlinear effects arising as a consequence of the procedure adopted in simulations S2-S6. The method used 
in simulations S7-S9 is clearly more suitable to compare the nonlinear effects of interacting DE models, and shows how the velocity-dependent acceleration 
(blue) is the most relevant mechanism in the suppression of small-scale power in these models. 



5.2 Baryon bias 

We now analyze the impact of the different individual effects of 
the DE-CDM interaction on a quantit y that has b een called in 
the literature the "gravitational bias" tealdi et al.K 2Q10). repre- 
sented by the ratio of the baryon to CDM power spectra defined as 
R{k, z) = Pbik, z) /Pc{k, z). This quantity shows how the density 
perturbations of baryons and CDM, which start with the same rela- 
tive amplitude in the early universe, evolve during structure forma- 
tion processes. Clearly, in absence of hydrodynamic forces acting 
on baryons (as it is the case for the numerical works under inves- 
tigation) the only difference in the growth of density perturbations 
between baryons and CDM is given by the DE-CDM interaction, 
such that for ACDM one expects R{k, z) = 1 \lk, z. This is indeed 
the case if one looks at Fig.|4l where the ratio R{k, z) is plotted 
at different redshifts for all the simulations under study. Again, up- 
per panels are for the simulations where the individual effects have 
been switched off right from the start, while the lower panels are 
for the case of switch-off only for z ^ z-ni. 

Also in this case the upper panels show a much larger scatter 
of the different simulations as compared to the lower panels, due to 
the large linear effects arising from an early suppression of the in- 
teracting DE effects. However, it is particularly interesting to notice 
in both types of approaches how the fifth-force term, that showed 
very little effects on the total power spectrum shown in Fig. [3] and 
discussed in the previous section, is now the leading mechanism 
in determining the evolution of the gravitational bias both at linear 
and nonlinear scales, such that its suppression significantly shifts 
up the value of R{k, z) towards the ACDM value of 1.0. Since the 



evolution of the gravitational bias is tightly related to the baryon 
fra ction of massive c ollapsed objects as galaxy clusters (as shown 
in iBaldi et al.ll2oTob . our analysis shows how the fifth-force term, 
far from being an irrelevant effect for interacting DE models, could 
actually determine observational features in the baryonic budget of 
groups and clusters of galaxies. 

It is then also very interesting to notice the impact of the 
velocity-dependent acceleration. Although the upper panels of 
Fig.m due to the large scatter induced by the high redshift suppres- 
sion of the interacting DE effects, do not allow a clear detection of 
this behavior, the lower panels distinctively show how the velocity- 
dependent acceleration differently affects the linear and the non- 
linear regimes of structure formation, as discussed in Sec. [3] In 
fact, while in the linear regime probed by the largest scales in the 
simulation the effect of the velocity-dependent acceleration is to en- 
hance the growth of CDM density perturbations, in the small scale 
regime it has the opposite effect of reducing nonlinear power. As a 
consequence, the suppression of this acceleration (that acts only on 
CDM particles) has the effect of increasing the value of the gravi- 
tational bias at large scales towards the ACDM value as compared 
to the full interacting DE model, and of further reducing it at small 
scales, as clearly represented by the crossover of the red line (full 
simulation) and the blue line (simulation S8) at fc ~ 1 h Mpc~^ 
in the last panel of Fig. |4] This behavior distinctly shows that the 
velocity-dependent acceleration cannot be treated in the same way 
in the linear and nonlinear regimes of structure formation, and that 
a suppression of this term already at high redshifts will inevitably 
determine a superposition of its linear and nonlinear impact onto 
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Figure 4. (Color online) The evolution of the gravitational bias R{k, z) as a function of inverse scale k for the ACDM model (black) and the interacting DE 
scenario (red). Upper plots show the results of the test simulations S2-S6 with high-redshift suppression of the individual effects, while lower panels show 
the results of simulations S7-S9 with suppression only for z ^ z^i = 2. The leading role of the fifth force (green) in determining a gravitational bias clearly 
appears in both sets of simulations. 



the final results of the simulations, making it very difficult to dis- 
entangle the two different contributions. 

5.3 Halo Mass Function 

We now compare the impact that each of the effects under study has 
on the halo mass function at 2; = and at z = 1. To do so, we first 
identify the halos formed within each of the simulations described 
in Table|4]by means of a Friends-of-Friends algorithm with linking 
length A — 0.2 x d, where d is the average particle spacing. We 
then i dentify halo substruc tures by means of the SUBFIND algo- 
rithm jSpringel et al.ll200ll) . and we compute the cumulative CDM 
halo mass function as a function of the halo virial mass Af200- The 
result is shown in Fig.|5]for the test simulations with a high-redshift 
switch-off of the different interacting DE effects (upper panels) and 
for the test simulations with Zni ~ 2 (lower panels). 

Also in this case, as expected, suppressing the effects of the 
DE-CDM interaction at high redshift determines a large scatter of 
the different test simulations at recent epochs, due to the signifi- 
cantly different linear growth history of the various runs. On the 
contrary, a low-redshift suppression of the individual effects pro- 
duces basically no scatter in the halo mass functions, which roughly 
follow the behavior of the full interacting DE model (red line) that 
by construction is normalized in order to give a comparable distri- 
bution of structures to the corresponding ACDM model (black line) 
at z = 0. 

The comparison of these two different scenarios is particularly 
relevant for the considerations that we will make in the next Sec- 
tion. In fact, it clearly appears from the upper panels of Fig.|5]that 
all the simulations with high-redshift suppression of the interacting 



DE effects will present very different structures at the present time, 
with significantly different masses and with possible large offsets 
in the final location of the individual halos, and this deviation from 
simulation to simulation is clearly more prominent at large masses. 
This will make particularly difficult to directly compare the proper- 
ties of individual halos within this set of simulations, since it will be 
more difficult to identify two halos formed in different simulations 
as the same object, as we will discuss in detail in the next Section. 

This again shows the potentially misleading superposition of 
linear and nonlinear effects that characterizes the approach of sup- 
pressing the new physical effects of interacting DE already at high 
redshifts in the simulations. The most suitable approach for study- 
ing the nonlinear re gime of interacting DE models is therefore the 
one first adopted bv lSaldi et alj 1 I2OICII) that allows to clearly disen- 
tangle the impact that each of the interacting DE effects has on the 
inner dynamics of nonlinear objects from the overall linear growth 
history. In fact, as we will show in the next Section, although the 
simulations displayed in the lower panels of Fig. |5] have roughly 
the same statistic distribution of massive halos, the internal proper- 
ties of individual halos still get significantly affected by each of the 
specific effects under investigation. 

5.4 Halo density profiles 

The last and probably most interesting test of our interacting DE 
models concerns the internal dynamics of highly nonlinear struc- 
tures as the massive CDM halos that form in each of the different 
numerical realizations under study. Since all the simulations are 
started from the same random realization of the linear power spec- 
trum in the initial conditions, it is expected that structures will form 



© 2010 RAS, MNRAS OOO.mfTil 



10 M. Baldi 





Figure 5. (Color online) The cumulative mass functions of CDM halos in ACDM (black) and interacting DE (red) models at z = 1 and 2 = 0. The large 
scatter of the different halo mass functions in the upper panels is due to the high-redshift suppression of the individual effects of the DE-CDM interaction in 
the simulations S2-S6. The same scatter does not appear if a low-redshift suppression is adopted. 



in the same position in all the simulations and this should allow a 
direct comparison of individual halos that formed in each of the 
runs. However, this correspondence cannot be exact due to the dif- 
ferent physical effects included in each simulation, and some offset 
in the final location of the halos is expected as a consequence of the 
specific dynamical evolution of the systems in each run. 

It is nevertheless very important to stress here that if mod- 
els are normalized in order to have the same linear perturbations 
amplitude at the present time, the offset in the position of the ha- 
los will generally be a fraction of their virial radius -R200 , and the 
local environment, the formation and merging history, and the dy- 
namical state of two corresponding halos in different s imulations 
will be in general very similar. This has been tested bv lBaldi et al.l 
where it has been shown that roughly 40% of the 200 most 
massive halos in a set of four hydrodynamical N-body simulations 
with different values of the DE-CDM coupling /? showed an offset 
smaller than one virial radius with respect to the location of their 
corresponding ACDM halo. 



On the contrary, if different models have a significantly differ- 
ent normalization of the linear perturbations amplitude at 2: = 0, 
the final structures formed in each simulation will significantly dif- 
fer from run to run, and will very likely be considerably offset from 
each other, as we showed in the previous section by comparing the 
mass functions of halos at 2; = and z = I'm the simulations with 
a high-redshift suppression of the different interacting DE effects 
(upper panels of Fig.HJ. 

Here we will apply a sim ilar identification criterion to the one 
adopted in iBaldi et al.l J201(]h . and within a given set of simula- 
tions we will identify with each other only those halos for which 
the location of the most bound particle (taken as the center of the 
halo) lies within a radius Ri — 7)i?200 from the center of a cor- 
responding ACDM structure. We can directly confirm our expec- 
tations by counting which fraction of halos can be identified with 
a given ACDM object by using different values of rj for the iden- 
tification radius Ri, first within the set of simulations that do not 
provide a common normalization of the linear perturbations ampli- 
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Figure 6. (Color online) Matter density profiles of two massive halos in our sample for the ACDM (black) and interacting DE (red) models. As in the other 
figures, upper panels refer to the procedure which adopts a high-redshift suppression of the individual effects of the DE-CDM interaction, while lower panels 
refer to the suppression only for z ^ z^i = 2. The lai'ge scatter present in the former case does not allow a clear comparison of the relative importance of each 
individual effect on the nonlinear density profiles due to a superposition of linear and nonlinear effects. On the contrary, the lower panels clearly show that the 
ve locity-dependent a cceleration (blue) is the most important effect in driving the reduction of the inner overdensity of halos, in full agreement with the results 
of lBaldi et all <2010l) . 



tude at 2 = (simulations S2-S6) and then within the simulations 
that as a result of a low redshift suppression of the individual inter- 
acting DE effects determine a substantially identical normalization 
at the present time (simulations S7-S9). 

For the former models, in fact, none of the halos at z = can 
be identified as being the same object in all the simulations when rj 
is set to unity. This fraction increases to 1% for r; = 2 and to 7% 
and 17% for = 3 and rj = 4, respectively. However, one should 
notice that an offset of 4i?200 corresponds for the most massive 
halos in our sample to a distance of 1.5 — 4.0 Mpc/h between the 
halos centers, which clearly does not guarantee that these halos can 
be safely considered as the same object forming in the different 
simulations. 

On the other hand, the latter models, which adopt the same 
procedure of Baldi et al. (2010), already have a fraction of 80% of 
identifiable halos for ri — 0.5, which increases to 91% for 77 = 1. 



It is therefore clear that these models will provide a much more 
homogeneous and statistically significant sample of halos for a di- 
rect investigation of the purely nonlinear effects of interacting DE 
models. 

This is confirmed by having a look at Fig.|6l where the density 
profiles of two halos in our sample are plotted for the S2-S6 sim- 
ulations with high redshift suppression of the specific interaction 
effects (in the upper panels) and for the S7-S9 simulations with 
suppression only for z ^ z^i- Evidently, the upper plots show a 
large scatter of the characteristic overdensity of the halos through- 
out the different realizations, and a significant discrepancy in the 
total mass, as it clearly appears from the large radii behavior of the 
density profiles in both plots. This confirms once again that for this 
type of simulations we are actually observing a superposition of 
linear and nonlinear effects even when considering the most highly 
nonlinear objects in our simulation box. It is then true, as noticed 
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bv lLi&Barro^ ( l2010bl) . that the suppression of the EDE compo- 
nent in the Hubble function (cyan line) and of the mass decrease of 
CDM particles (dotted orange line) determine the most significant 
increase (even with respect to the velocity-dependent acceleration) 
of the halo overdensity in the inner regions, but the same increase 
is visible at all radii, as it clearly appears when looking at the be- 
havior of these models in the outskirts of the halos. This however 
does not happen for the case where the velocity-dependent term is 
switched off (blue line), which shows a comparable characteristic 
overdensity at the virial radius to the ACDM model (black line) 
and to the full interacting DE simulation (red line). This superpo- 
sition of linear and nonlinear effects can then lead to the wrong 
co nclusion - claimed as a partial discrepancy with the early results 
of lBaldi et alH2010l) - that the background evolution and the mass 
variation of CDM particles are the most important features of inter- 
acting DE models in the nonlinear regime. 

This is clearly not true if we now consider the two bottom 
panels of Fig. [6] where the same halos are compared for the sim- 
ulations where the individual effects are swhitched off only at low 
redshift. As we already stressed above, this produces halos with a 
much more similar environment and with a comparable formation 
redshift, and allows a much clearer identification of corresponding 
halos in the simulation set. In fact, from the bottom panels of Fig.|6] 
one can notice that the halos in all the simulations have the same 
characteristic overdensity at the virial radius, nevertheless show- 
ing very interesting departures from the ACDM and from the full 
interacing DE simulation in the inner regions. Therefore, this nor- 
malization allows to test the nonlinear behavior in the halo core 
in an almost independent way from the linear normalization of the 
surrounding environment. 

Using this more meaningful procedure for the comparison 
of the different effects of the DE-CDM interaction s hows r esults 
in full agreement with what found by iBaldi et al.l feOld) . The 
velocity-dependent acceleration (blue line) is clearly the most im- 
portant effect in lowering the inner overdensity of massive halos, 
such that its suppression determines the larg est overdensity in - 
crease in the halo cores. As already found by iBaldi et al] ( I2OI0I) . 
the mass variation has a minor but still significant effect, while the 
fifth force shows practically no impact on the internal dynamics of 
the halos. 

These results therefore fully address the apparent discrepancy 
claimed bv lLi&Barr"owl ( l2010bl) that found evidence for the oppo- 
site hierarchy of the velocity-dependent acceleration and mass vari- 
ation effects in lowering halos inner overdensities. The discrepancy 
only arises as a consequence of the n umerical procedures adopted 
in the study of iLi & Barrowi l l2010bh that determines a superposi- 
tion of linear and nonlinear effects of the velocity-dependent accel- 
eration, and a superposition of "mass constancy" (i.e. nic = 0) and 
"mass normalization" (i.e. the effective Sic) for what concerns the 
variable mass of CDM particles. 



6 CONCLUSIONS 

Interacting DE models have become very popular in recent years 
as a viable alternative to the standard ACDM cosmological model. 
The investigation of these alternative cosmologies has shown sig- 
nificant progresses by moving from the study of the background 
evolution to the analysis of linear perturbations effects and ulti- 
mately to the impact of interacting DE models on nonlinear struc- 
ture formation by means of specific modifications of N-body algo- 
rithms. 



The main effects through which an interaction between DE 
and CDM can affect the growth of density perturbations range from 
a modified background expansion, to a time variation of the CDM 
particle mass, to a long-range attractive "fifth force" between CDM 
particles and a "modified inertia" in the form of a new velocity- 
dependent acceleration. 

In this work we have presented a detailed study of the impact 
that interacting DE models have on linear and nonlinear structure 
formation processes. Our work significantly extends previous anal- 
yses and aims at a direct comparison between different numerical 
procedures previously adopted in the literature to investigate the 
relative importance of each of the above mentioned physical effects 
that characterize interacting DE cosmologies. 

By means of a suitable modification of the N-body code 
GADGET-2, we have performed a series of collisionless cosmologi- 
cal N-body simulations for a standard ACDM model and for an in- 
teracting DE model with coupling /3 = 0.24. Such a large value for 
the coupling is already ruled out by several observational probes, 
however our aim here was not to describe an observationally vi- 
able scenario but rather to explore in detail how an interaction be- 
tween DE and CDM affects structure formation, and a large cou- 
pling clearly makes this task more easily achievable by amplifying 
the effects under investigation. 

In addition to the ACDM and the interacting DE simulations, 
we have carried out other 8 simulations in which each of the spe- 
cific effects of the DE-CDM interaction has been in turn artificially 
suppressed, in order to quantify its relative contribution to the 
different peculiar features of interacting DE cosmologies. In doing 
so, we have applied two different procedures previously adopted in 
the literature, namely the selective suppression of each individual 
effect only at the latest stages of structure formation w hen most of 
the no nlinear processes take place (as first proposed in lBaldi et al.l 
l2010l) . and the switch off of the different ef fects during the whol e 
simulations (as more recently done by iLi & Barrowi l2010bh . 
Our analysis therefore allows a direct comparison of these two 
methods and provides a direct way to test the relative importance 
of the different features of interacting DE models at high and low 
redshifts. 

As a first test we have studied the impact of each individual 
effect of the DE-CDM interaction on the evolution of the matter 
power spectrum at different redshifts. The global effect of interact- 
ing DE on the matter power spectrum is to suppress power at small 
scales while the large scale amplitude has the same normalization 
as for ACDM at 2; = 0. This first test already shows a signifi- 
cant difference between the two numerical procedures mentioned 
above. In fact, the suppression of any of the effects of the inter- 
action right from the start of the simulations determines a differ- 
ent linear growth of density perturbations that produces a signifi- 
cant scatter in the large scale normalization of the power spectra 
at 2 = 0. As a consequence, the final power spectrum of each of 
these simulations is modified at small scales by a superposition of 
the linear and nonlinear impact of the suppressed effect, making a 
direct test of the relative importance of each effect quite difficult. 

On the contrary, if the suppression is limited to low redshifts, 
the large scale normalization of the power spectrum remains 
consistent with the original ACDM and interacting DE models, 
and the small scale differences from simulation to simulation 
will be due only to the nonlinear impact of each specific effect. 
With this more suitable comparison procedure we have clearly 
shown that the velocity-dependent acceleration of CDM particles 
is the most important effect in suppressing small scale power. 
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with the mass variation playing a minor but still significant role, 
while the fifth force has basically no effect in the nonlinear regime. 

We have also tested how each of these effects alter the relative 
evolution of density perturbations in the uncoupled baryonic 
component as compared to CDM. In this case, the hierarchy of the 
different effects is reversed with respect to the impact on the power 
spectrum, and the fifth force shows the most important contribution 
to the faster growth of CDM perturbations with respect to baryonic 
perturbations both at linear and nonlinear scales. This different 
growth gives rise to the so called "gravitational bias" between the 
two components which determines a significant baryon depletion 
of collapsed objects as compared to ACDM. Therefore, the fifth 
force is found to be the driving mechanism for the reduced 
baryon fraction of massive halos that characterize interacting DE 
cosmologies. 

We have then investigated the impact on the halo mass func- 
tion of each of the effects under study, where once again the two 
different numerical procedures used in our analysis show signifi- 
cant differences. Also in this case, suppressing individual effects 
from the beginning of the simulation determines a large scatter in 
the statistical distribution of collapsed objects at z = with respect 
to the original ACDM and full interacting DE simulations that by 
construction give very similar mass functions at the present time. 
This scatter is maximum at large masses and shows how, in the sim- 
ulations run with this procedure, structures cannot be expected to 
form in the same locations and have similar masses and formation 
histories from simulation to simulation. This clearly makes any at- 
tempt to directly compare individual halos in different simulations 
particularly hard, since it will be difficult to safely identify objects 
in different runs as being the same structure, as we have directly 
shown in our analysis. 

On the contrary, again, we showed that suppressing individual 
effects of the DE-CDM interaction only at the latest stages of 
structure formation produces very similar statistical distributions 
of halos at the present time, with basically no scatter around the 
original mass functions of ACDM and of the full interacting DE 
model. This procedure therefore results clearly more suitable also 
in order to perform direct comparisons between individual halos in 
different numerical realizations, as the spatial locations of bound 
objects as well as their local environment and formation histories 
will be very similar from simulation to simulation. 

We have finally studied the relative importance of each indi- 
vidual effect of the DE-CDM interaction on the radial density pro- 
files of CDM halos. Consistently with the results found for the halo 
mass functions, we have shown that a suppression of any specific 
effect from high redshifts determines a large scatter of the chara- 
cteristic overdensity of the halos at z = 0, which have a very dif- 
ferent total mass and show significant differences in the amplitude 
of the density profiles at all radii. This behavior again witnesses the 
superposition of linear and nonlinear effects when this numerical 
procedure is adopted. 

On the other hand, switching off the individual effects only 
at low redshifts preserves the characteristic overdensity of corre- 
sponding halos that show the same amplitude of their density pro- 
files at large radii. As a consequence, the impact of each individual 
effect on the internal dynamics and on the distortion of the profile 
at small radii is much more clearly visible and can be safely disen- 
tangled from the linear normalization of the environment in which 
each halo is embedded. By adopting this more suitable procedure. 



we found once again that the most relevant effect in reducing the 
inner overdensity of halos in interacting DE models is given by 
the velocity-dependent acceleration of CDM particles, whereas the 
mass variation has a minor but not negligible impact, while the fifth 
force shows no influence whatsoever in this highly nonlinear con- 
text. 

These results are in full agreement with what previously found 
and show that the discrepancies claimed by 



bvlBaldi etid 
iLi & Barrowl I 



201 Oh) are essentially due to the different procedure 



adopted in their numerical study, which determines a superposition 
of linear and nonlinear effects and is therefore less suitable to 
compare interacting DE models in the highly nonlinear regime of 
structure formation. 

To conclude, we have performed a wide and detailed study of 
interacting DE cosmologies aimed at comparing the relative impor- 
tance that each of the new physical effects that characterize these 
models has in affecting linear and nonlinear structure formation 
processes. Our study significantly extends previous works and al- 
lows a direct comparison between different numerical procedures 
recently adopted in t he literature . Our outcomes fully confirm the 
early results of Baldi et alj j201Q). and sho w that the apparent dis- 
crepancies found bv lLi & Barrowl bOlObh are actually due to the 
different numerical setup adopted by the latter work, which in- 
evitably determines a superposition of linear and nonlinear effects 
and is therefore not particularly suitable to compare the different 
characteristic features of interacting DE models in the nonlinear 
regime. 
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